Energy distribution and cooling of a single atom 
in an optical tweezer 
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We investigate experimentally the energy distribution of a single 87 Rb atom trapped in a strongly 
focused dipole trap under various cooling regimes. Using two different methods to measure the mean 
energy of the atom, we show that the energy distribution of the radiatively cooled atom is close 
to thermal. We then demonstrate how to reduce the energy of the single atom, first by adiabatic 
cooling, and then by truncating the Boltzmann distribution of the single atom. This provides a 
non-deterministic way to prepare atoms at low microKelvin temperatures, close to the ground state 
of the trapping potential. 

PACS numbers: 37.10.Jk, 42.50.Ct, 03.67.-a 



I. INTRODUCTION 

In the last thirty years, the manipulation of single 
quantum objects has been the subject of considerable at- 
tention. After the pioneering experiments of H. Dehmelt 
with a single electron in a Penning trap a single ion 
was trapped in a Paul trap 0, which led to the obser- 
vation of quantum jumps [31,14]. Several techniques were 
then developed to cool the single trapped particle, such 
as optical sideband cooling, which was experimentally 
demonstrated in 1989 ||. More recently, the interest in 
the manipulation of single quantum objects has grown 
with the recognition that they can carry quantum infor- 
mation, and that they therefore can be a resource for 
quantum information processing [||. In this context, sin- 
gle trapped laser cooled ions have accumulated an im- 
pressive amount of results in recent years (see e.g. Q). 

The trapping of single neutral atoms, on the other 
hand, is more recent, owing in part to their weaker inter- 
action with the electromagnetic field. So far, the trap- 
ping of neutral atoms has been realized in tig ht magneto- 
optical traps [1, [§] or optical dipole traps [l(| [HI, ED, • 
As the trapping potentials are typically less than a mil- 
liKelvin deep, laser cooling is the usual starting point of 
the manipulation of single atoms. As is the case for ions, 
it is tempting to further reduce the energy of the trapped 
atom, to ultimately reach the ground vibrational state 
of the trapping potential. In the framework of quan- 
tum computing, for example, the entanglement of two 
atoms via controlled collisions usually requires ground 
state cooling (see e.g. [Hj]). 

The knowledge and control of the energy of a single 
atom, trapped in a sub-micron optical dipole trap, is im- 
portant for many purposes. For instance, in a recent ex- 
periment [l5j . we have used two single rubidium atoms, 
trapped in two optical tweezers separated by a few mi- 
crons, as single-photon sources. When the two photons 
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(each photon coming from one atom) are superimposed 
on a beamsplitter, we observe a "coalescence" effect : 
the two (indistinguishable) photons copropagate from the 
same output of the beamsplitter, due to a two-photon 
interference effect. We have shown that the ultimate vis- 
ibility of this effect is limited by the spread in frequencies 
of the emitted photons, which in turn is due to the spread 
in energy of the emitting trapped atoms. In another ex- 
periment [l6| , using Raman transitions we have prepared 
a single atom in a superposition of two internal states of 
the hyperfine manifold, and we have studied the dephas- 
ing of this quantum bit. In this case, due to the fact that 
successive atoms have different energies, each realization 
of the qubit evolves differently with respect to a local os- 
cillator. This results in a loss of coherence, that we were 
nevertheless able to compensate for by applying rephas- 
ing "spin-echo" techniques. For this application, a lower 
energy spread of the trapped atom would also increase 
the dephasing time, thus avoiding the extra "spin-echo" 
sequence. 

In the present paper, we analyze in more detail the 
methods used in these previous works to determine the 
mean energy of a single atom in the optical tweezer 
[TEl [ilil . fl?! . Similar or complementary techniques have 
been reported by other groups [12|, [l8| . We then exploit 
these temperature measurements to characterize several 
cooling methods to further decrease the temperature of 
the atoms, and we show how to approach the ground 
state level in the two radial directions of the trap, which 
exhibit the strongest confinement. 

The paper is organized as follows. We describe in sec- 
tion |TT] the experimental setup. In section IITTI we present 
a "release and recapture" technique used to measure the 
temperature of the atoms. In section IIVI we develop 
a cooling sequence used to laser-cool the atom in the 
tweezer. In section IVl we perform the spectroscopy of the 
energy distribution of the atom in the tweezer and find 
that this distribution is very close to a Boltzmann distri- 
bution. Section IVT1 describes how the mean energy of the 
single atom can be decreased by truncating this Boltz- 
mann distribution. This provides a non-deterministic 
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FIG. 1: (Color online) Schematic of the experimental setup. 
A beam of 87 Rb atoms, originating from an oven, is slowed 
using a Zeeman coil (not shown) and counter propagating 
slowing (fii) and repumping (f2s) lasers at ~ 780 nm. An 
optical molasses is formed at the intersection of six counter 
propagating lasers (£^3) and repumpers (£^4) operating also 
at ~ 780 nm. An aspherical lens is used to tightly focus 
a ~ 850 nm laser (fis) in the region of the optical molasses. 
The fluorescence of the atom, due to the molasses and Zeeman 
lasers, is collected using the same aspherical lens. A dichroic 
mirror together with other spatial and spectral niters are used 
to filter the fluorescence signal, which is detected using single- 
photon counting avalanche photodiode (APD). 



method to prepare a single atom close to the ground state 
of the trapping potential. In section IVIII we reduce fur- 
ther the temperature of the atoms by adiabatically lower- 
ing the potential trap. Finally, we conclude by discussing 
some implications of these results. 

II. EXPERIMENTAL SETUP 

Figure [T] shows a schematic of the experimental setup. 
A beam of 87 Rb atoms is slowed using a Zeeman slowing 
technique (l9j |. The counter- propagating laser beam at 
~ 780 nm (denoted f2i) is 18.7 T red detuned from the 
F = 2 — > F' = 3 cycling transition, and combined with 
a repumping laser (fig) that is 3.7 T red detuned from 
the F = 1 — > F' — 2 repumping transition. Here, T « 
2ir x 6 x 10 6 rad/s is the linewidth of the D2 transition. 

We form an optical molasses at the intersection of six 
counter-propagating cooling lasers that are 4.2 T red de- 
tuned from the F = 2 — ► F' = 3 transition (denoted f^), 
together with six repumping lasers (f2 4 ) resonantly tuned 
to the F = 1 — > F' = 2 transition. This optical molasses 
is used as a reservoir of cold atoms from which the dipole 
trap is loaded. 

We produce the optical dipole trap by focusing a 
~ 850 nm laser in the center of the molasses, using a 
large numerical aperture (NA = 0.5) aspherical lens p7| . 
The optical dipole trap has a measured optical waist of 
w = 1.03 ± 0.01 (im. A comprehensive characterization 



of the optical system and of the dipole trap itself is pre- 
sented in detail elsewhere [l7|. Assuming a Gaussian 
beam profile in the transverse direction for the dipole 
trap and a Lorentzian profile along the propagation axis, 
we calculate the radial and axial frequencies and find 
u± ~ 160 kHz and v\\ ~ 30 kHz respectively for 10 mW 
of laser power. The potential trap depth, proportional 
to the laser power, is U ~ 2.8 mK for a 10 mW of laser 
power. 

This optical tweezer allows us to trap single Rb 87 
atoms via a collisional blockade mechanism, which pre- 
vents two or more atoms from being trapped simulta- 
neously due to inelastic collisions [lj, [l7], H(J. The 
~ 780 nm fluorescence of the trapped single atoms is 
detected using the same aspherical lens that is used to 
create the dipole trap. This fluorescence is filtered using 
a dichroic mirror, together with other spatial and spectral 
filters, and then detected using a single photon counter 
(Avalanche Photodiode), with an overall collection and 
detection efficiency that is estimated to be ~ 0.5%. From 
the output signal of the single-photon counter we distin- 
guish between the presence (or absence) of an atom in the 
tweezer, within a 10 ms acquisition period, with a confi- 
dence better than 99%. Each experimental sequence that 
will be presented in this paper is triggered on a detected 
single atom in the dipole trap. 

Our goal in the following sections will be to investigate 
the energy distribution and temperature of the atoms 
trapped in the optical tweezer, under various cooling con- 
ditions. 



III. RELEASE AND RECAPTURE TECHNIQUE 

In this section we will describe a "release and recap- 
ture" method that we use to determine the temperature 
of the atoms in the dipole trap [HI, [l(| • We emphasize 
that what we mean by the temperature of the atoms is 
a temperature extracted by averaging the energy over 
many realizations of the same experiment with a single 
atom. The premise of this method is that information 
can be obtained about the energy distribution of a single 
atom by switching off the dipole trap for a variable time 
At (as illustrated in Fig.[^a)) and then determining the 
probability of recapturing the atom, which is denoted as 
P rr (Af). On average, an atom with a high energy is more 
likely to escape during the release time At, compared to 
an atom with a low energy, as is well-known in the con- 
text of laser cooling of atomic samples [2l|, ■ In the 
single atom regime used here, this method is particularly 
robust, since on each repetition of the experimental se- 
quence the result is binary : either the atom has escaped 
the trapping region, or it is recaptured. This binariza- 
tion of the measurement outcome makes it immune to 
various noise sources, in particular to fluctuations in the 
fluorescence level of the atom, which is used to decide on 
the presence of the atom. 

Figure (Ha) illustrates a schematic of the "release and 
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FIG. 2: (Color online) Measurement of the temperature of 
the atoms using the "release and recapture" method, (a) 
Schematic of the experimental sequence (see text), (b) Ex- 
perimental results showing the probability of recapturing the 
atom, P r r(At), as a function of release time At. (i) (Circles) 
Temperature of the atoms just after the loading in the dipole 
trap. Each data point is the accumulation of 100 sequences. 
Superimposed on this data is a fit by the Monte-Carlo sim- 
ulation of the "release and recapture" method, which is the 
average of 200 trajectories for each release time. The tem- 
perature is 168 ± 6 /jK for a trap depth of ~ 2.8 mK. (ii) 
(Diamonds) For comparison we have also plotted the results 
for the temperature measurement using the same method, 
but after the atom has been laser-cooled in the tweezer. Each 
data point is the accumulation of 200 sequences. The best-fit 
temperature is 31 ± 1 /iK for a trap depth of ~ 2.5 mK. 



recapture" sequence, which goes as follows, (i) After the 
atom is trapped, the cooling and loading lasers (fix, O2, 
0$ and O4) are switched off. (ii) The dipole trap laser 
(O5) is then turned off for a time At. (iii) Finally, we per- 
form a fluorescence detection by switching all the cooling 
and loading lasers back on, thus determining whether the 
atom has escaped or not. Figure [21 shows typical experi- 
mental results for the "release and recapture" technique, 
where the probability of recapturing the atom is plot- 
ted as a function of the release time At (more details 
about this experiment will be presented in the following 
section). 

In order to extract the temperature from these mea- 
surements, we perform Monte-Carlo simulations of the 
trajectories of the single atoms, to determine the proba- 
bility that the atoms are recaptured after a release time 
At. This simulation requires the knowledge of the en- 
ergy distribution of the atom in the trap. We assume 
that we are in the harmonic approximation of the trap- 
ping potential, and that the position- velocity distribution 
of the atom follows a thermal Maxwcll-Boltzmann law. 
The standard deviations of the position of the atom in 
the axial and radial directions of the trapping potential 
are Ax\\ — k^T / 'muj 2 and Ax± = \J k&T / moj' 2 ^ re- 
spectively, while th e standard deviation of the velocity 
is Av — y k^T '/m, where T is the temperature. The 
parameters of the optical trapping potential have been 



determined previously [lTf , and take into account the tilt 
due to gravity, which creates a local maximum in the po- 
tential along the gravity axis. The effective trap depth is 
defined as the height of this potential barrier with respect 
to the minimum of the potential. 

The simulation goes as follows: using the 
Maxwcll-Boltzmann distribution we randomly gen- 
erate the three-dimensional position-velocity vector 
(xi, yi, Z{, i> x .i, Wy.i, w z .i) of the atom in the potential for a 
given temperature. We then calculate the trajectory of 
each atom during the release time At in a single time-step 
calculation. The position of the atom after this release 
time is {x\ + v X}i At,yi + v y _[At — gAt 2 /2,Zi + v Z: iAt), 
where the potential gradient due to gravity along the 
y-axis has been included. We finally determine the 
energy of each atom when the trap is turned back on. If 
the total energy of the atom after the time-of-flight is 
smaller than the effective trap depth, we consider that 
the atom is recaptured, otherwise it is considered to 
have escaped. 

By simulating many trajectories of the single atoms for 
each release time At, we can numerically determine the 
probability of recapturing the atom as a function of this 
release time. These simulations are repeated for a range 
of different temperatures. We calculate the weighted 
least square value (% 2 ) between the experimental data 
and the simulated results for each temperature. The chi- 
square value is defined as \ 2 — — Pn(Ati)] 2 /of} 

where P lr (At) is the experimental data that has an un- 
certainty a and / is the value predicted by the model. 
The best-fit temperature corresponds to the Boltzmann 
distribution that minimizes this least square value. 

We performed the "release and recapture" technique to 
determine the temperature of the atoms just after they 
were loaded in the trap. Figure [U^b) shows the Monte- 
Carlo simulation that best fits the corresponding exper- 
imental data. In this example, the best fit temperature 
is T = 168 ± 6 in a trap depth of U ~ 2.8 mK. The 
ratio between the trap depth and the temperature of the 
atoms is 15, indicating that we are approximately in the 
harmonic regime of the trapping potential. 

Due to experimental miscounts and losses, the prob- 
ability of recapturing an atom for zero release time was 
less than unity, being instead typically 0.95. We scaled 
the Monte-Carlo simulation by this factor to achieve a 
better fit to the experimental results. The uncertainty 
on the temperature that minimized the chi-square value 
is defined as ctt = \/2(d 2 x 2 /<9T 2 ) -1 , assuming that the 
X 2 function varies harmonically around this minimum 
point [23]. Furthermore, we have included in the error 
budget for the temperature, the statistical uncertainty in 
fitting a parabola function to this chi-square minimum. 

We will use this "release and recapture" technique 
throughout the rest of the paper as a diagnostic to mea- 
sure the temperature of the atoms in the optical tweezer. 
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FIG. 3: (Color online) Schematic of the temporal sequence 
for the laser cooling sequence. The details are given in the 
body of the text. The box labeled Temperature Experiments 
shows the point in the sequence where we can manipulate 
the intensity of the dipole trap laser in order to measure the 
temperature of the atoms. (See Sections IVl IVII and IVII|I . 



mental results after the laser-cooling, together with the 
best- fit simulation results. This corresponds to a temper- 
ature of 31 ± I /jK for a trap depth of ~ 2.5 mK, which 
yields the ratio U/T ~ 81. This represents more than 
a factor 5 reduction in the temperature with respect to 
the temperature of the atoms directly loaded from the 
optical molasses. 

We note that the optimized cooling sequence is 
achieved when the laser detuning is brought closer to 
resonance with respect to the atomic transition. This 
fact may seem inconsistent with the sub-Doppler theory 
of laser cooling (see e.g. (HI), which predicts that the 
temperature is lower for larger detuning. However one 
has to take into account that when the atom is trapped 
in the tweezer, the resonance frequency is light shifted 
by about 10T with respect to the free space case. There- 
fore the detuning experienced by the trapped atom varies 
only from — 14.5T to — 12I\ The extra cooling effect is 
then attributed to the decrease (by a factor 3) in the in- 
tensity of the laser beams, rather than to the frequency 
change, now in reasonable agreement with sub-Doppler 
laser cooling. 

This optimized cooling technique is the starting point 
for each experiment presented throughout the rest of the 
paper. We will now turn our attention to validating the 
assumption about the thermal distribution. In addition, 
we would like to cross-check the temperature results using 
an alternative method. 



IV. RADIATIVE COOLING OF A SINGLE 
ATOM 

In this Section, we discuss the laser cooling of a sin- 
gle atom after it has been captured in the dipole trap. 
For that purpose, the cooling sequence was optimized by 
controlling the molasses lasers parameters (cooling time, 
frequency detuning, intensity), as well as the trapping 
laser intensity. This was done by maximizing the recap- 
ture probability of the single atom after a "release and 
recapture" experiment with At = 10 fis. 

Figure [3] shows the optimized laser cooling sequence. 
Initially, all lasers are on and the molasses cooling lasers 
are detuned by Ai — — 4.5T to maximize the loading rate 
of the molasses and optical dipole trap, (i) The sequence 
is triggered on the detection of an atom, where the dipole 
trap laser power is ~ 10 mW. (ii) The Zeeman laser and 
the Zeeman repumper laser are switched off. At the same 
time, the cooling lasers are linearly detuned from Ai = 
— 4.5r to A2 = — 2T in ~ 1.2 ms. During this time, 
the intensities of the molasses lasers are reduced from 
Ix ~ 2 mW/beam, to approximately one third of their 
initial intensity. After this cooling period the cooling and 
repumping lasers are switched off. 

Using this optimized laser-cooling sequence we were 
able to significantly reduce the temperature of the atoms. 
Figure [2Kb) shows the "release and recapture" experi- 



V. ENERGY DISTRIBUTION OF A SINGLE 
ATOM IN AN OPTICAL TWEEZER 

In this section we determine the energy distribution of 
the single atom in the dipole trap. Provided that this 
distribution is indeed thermal, the corresponding tem- 
perature can be evaluated. To achieve this, we apply a 
type of spectroscopy method proposed and experimen- 
tally demonstrated by Alt et al. in ref. [l8l |. 

This method involves the adiabatic lowering of the trap 
depth to a point where the atom can potentially escape. 
As shown in Fig. [4ja), by lowering the trap depth, an 
atom with initial energy ~E\ in a trap of depth U\ will 
eventually have an energy E csc that is equal to the final 
trap depth U csc and hence the atom will escape the trap. 
Thus, measuring the depth at which the atom escapes 
the trap yields information about the initial energy of 
the atom. 

Figure H^b) shows a schematic of the experimental 
sequence for this method, which goes as follows: (i) 
The single atom is trapped and laser cooled, (ii) The 
dipole trap depth is then adiabatically lowered from 
Ui ~ 2.8 mK to t/ m in- The lowering time is approximately 
~ 2.5 ms. The trap depth is held constant for a dura- 
tion of ~ 20 ms, to allow an atom with an ener gy g reater 
than the trap depth sufficient time to escape [301 ]. (iii) 
After this waiting time, the trap depth is adiabatically 
raised to the initial trap depth U in ~ 2.5 ms. Finally, 
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the laser beams are switched back on to determine if the 
atom has escaped or not. We leave to section IVIII the 
discussion about the adiabaticity of this process. Here 
we assume that the hypothesis for the adiabatic change 
of the potential is fulfilled. Figure [4] (c) shows the ex- 
perimental probability of recapture for various minimum 
trap depths t/ m i n normalized to the initial trap depth XJ\ 

It is now necessary to map the minimum trap depth 
U m - m back to the corresponding energy in the initial 
trap Ei. To perform this mapping we calculate the one- 
dimensional action that can be expressed as the integral 
S(E, U) = J Xmax v / 2m[E -U(x)]dx, where E is the en- 
ergy, U(x) is the potential and x max is the position where 
the atom has zero kinetic energy. In the adiabatic limit, 
the action of the atom in the trap is conserved as the 
trap depth is adiabatically lowered [2{|. We solve the 
constant action equation S(Ei, U) = S(U CSC , U csc ) in the 
radial direction of our potential trap that includes the 
gravity potential gradient. The numerical results of the 
Ucsc/Ui — ► Ei/Ui mapping is plotted in Figure [4] (e). We 
subsequently apply the mapping U csc /Ui — ► E\j\J\ to our 
experimental data. Figure 0] (d) shows the percentage 
of atoms that have an energy less than £\ in the initial 
trap of depth Ui plotted as a function of the normalized 
energy Ei/Ui. 

We now wish to compare the measured energy distri- 
bution to a thermal one. In the harmonic trap limit (i.e. 
when the ratio of the trap depth to the temperature of the 
atoms satisfies Ui/k^T ^> 1), we calculate the Boltzmann 
energy distribution for a temperature T. We assume that 
the normalized energy distribution of the atom is 



ME) = 



1 



2(k B T) 



:E 2 e -E/k B T 



(1) 



where E 2 corresponds to the three dimensional density 
of states. The probability that an atom has an energy 
less than E is defined by P s , 
is given by 



,{E) = /o / th (£7') dE' and 



P SUIV (E) = 1 



1 



(2) 



where T) = E/k^T. Equation ((2|) can be interpreted as 
the survival probability of the atom remaining in the 
dipole trap after the truncation of the Boltzmann dis- 
tribution. 

Figure Ud) shows a fit of the integrated Boltzmann 
distribution P (Ei/Ui) to the experimental data using 
Eq. ([2|) . To achieve a better fit to the data, the maximum 
survival probability is rescaled to 0.95, to account for 
non-ideal recapture probability of our experiment when 
U m i n = U. We find that the measured energy distri- 
bution is well fitted by a thermal distribution. This 
fact shows that over many realizations of the experi- 
ments, the energy of successive atoms after laser cooling 
is well described by a Boltzmann distribution, in agree- 
ment with several previous theoretical and experimental 
studies I29| . 
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FIG. 4: (Color online) Spectroscopy of the energy distribution 
of a single atom using the adiabatic lowering method, (a) 
Schematic showing that due to the adiabatic lowering of the 
dipole trap, an atom with initial energy Ei in a trap depth Ui 
eventually has an energy equal to the trap depth, i.e. -Emin = 
Umin = Uesc- (b) Schematic of the experimental sequence (see 
text), (c) Experimental results of the survival probability for 
various minimum trap depths (7 m m. Each data point is the 
accumulation of 200 sequences, with error bars due to the 
binomial statistics. Solid line is a theoretical fit to the data, 
(d) Reconstruction of the Boltzmann distribution in the initial 
trap of depth Ui. Plot showing the percentage of atoms with 
an energy less than Ei (in the initial trap of depth Ui) versus 
the normalized energy Ei/Ui. Solid line is a theoretical fit to 
the data (see text). Dashed line represents the corresponding 
Boltzmann distribution (e) Numerical results for solving the 
constant action equation S(Ei, Ui) = S([4ac, U esc ). 
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The Boltzmann distribution in Figure H|d) has a cor- 
responding temperature of T = 33 ± 2 /iK for a trap 
depth Ui ~ 2.8 mK, where the error in the temperature 
corresponds the statistical uncertainty in the fit. This 
temperature is in good agreement with the "release and 
recapture" temperature results presented in the previous 
section. 

We point out that the adiabatic lowering of the dipolc 
trap depth is a powerful tool in manipulating the energy 
distribution of the atom. For example, at very low trap 
depths the more energetic atoms will escape with a higher 
probability. This adiabatic lowering indeed acts as a filter 
such that only the coolest atoms will remain in the trap, 
as we will now show. 



VI. REDUCTION OF MEAN ENERGY BY 
TRUNCATING THE BOLTZMANN 
DISTRIBUTION 

We have shown in Sections IIVI and [V] that we can 
cool the atoms down to a temperature of approximately 
T = 33 fiK using laser cooling techniques. It is possible 
to further reduce the mean energy of the trapped atom 
by truncating the Boltzmann distribution at an energy 
ftrunc- Here, we perform this truncation by adiabatically 
lowering the trapping potential, to a point where atoms 
with a higher energy are likely to escape, while atoms 
with a lower energy remain trapped. This filtering via 
the lowering of the trap potential appears to be similar 
to evaporative cooling, but a crucial difference is that 
there is no collision-induced rethermalization, since the 
experiment is done at the single atom level [27j |. 

Figure E[a) shows a schematic of the experimental se- 
quence that we use to investigate the energy distribution 
of the truncated Boltzmann distribution. This experi- 
mental sequence is identical to the one described in sec- 
tion [Vj but we now perform a "release and recapture" 
sequence directly after it. The entire temporal sequence 
is shown in the inset of figure E^b) . The additional "re- 
lease and recapture" sequence allows us to determine the 
temperature of the atoms after the truncation process 
and after the trap depth has been returned to the initial 
depth. 

To determine the truncation energy fTtrunc m the trap 
of potential depth Uf = U, we first calculate the trunca- 
tion energy in the shallow trap of potential depth U m - m 
(using the measured minimum optical dipole trap power) . 
We then use the numerical results that map t7 m i n — ► E\ 
(see Fig. HJe)) to determine the truncation energy J7 tr unc 
in the trap of potential depth Uf. 

In the case of a truncated Boltzmann distribution, it 
is not sufficient to assume that the position- velocity vec- 
tor of the atom follows a Gaussian Maxwell Boltzmann 
distribution, as we did for the "release and recapture" 
experiments in Sections IIIII and IIVI To account for this 
effect, in the following simulations we consider that the 
energies of the atoms are drawn from a truncated Boltz- 
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FIG. 5: (Color online) Investigation of a truncated Boltz- 
mann distribution using the "release and recapture" method, 
(a) Schematic of the experimental sequence, (b) Probabil- 
ity of recapture for zero release time P TI (At — 0) plotted as 
a function of the minimum trap depth U m i n . P TI (At — 0) 
is defined as the mean of the first four recapture probabili- 
ties with the smallest release times of a given data set (e.g. 
see Fig. EJb)), where the error bars are due to normal statis- 
tics. Solid line is a theoretical fit to the data. Inset shows 
a schematic of the temporal experiment sequence, (c) The 
temperature T plotted as a function of Umin- Solid line 
corresponds to a constant temperature of 33 fiK. Grey bar 
corresponds to the standard deviation of the spread in the 
measured temperatures for the minimum trap depth U m in ~ 
2.8 mK. (d) The calculated mean energy (E) of the trun- 
cated Boltzmann distribution shown as a function of U m i u . 
The solid line is a theoretical prediction of the mean energy 
of a truncated Boltzmann distribution with a temperature of 
33 fjK. 
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maim distribution corresponding to a temperature T and 
a truncation energy C/ tr unc- The Boltzmann distribution 
is discretized into N bins with energies £j equally spaced 
by AU = U tlunc /N, for AU/2 < E, < ([/ trunc - AC//2). 
For each distribution, the number of simulations with 
atoms with an energy Ej is weighted according to this 
discretized Boltzmann distribution. For an atom with 
energy Ej, we randomly distribute the energies among 
the three directions given by the cartesian axes. The pe- 
riod of motion of the atom is then randomly chosen for 
each of the three cartesian axes (ensuring that the total 
energy of the atom remains Ej). Finally, the trajectories 
of the atoms are simulated during the time of flight as de- 
scribed in Section IIIII We typically use several hundred 
atom trajectories for each release time At, and we dis- 
cretize the Boltzmann distribution into 10 bins. When 
the truncation energy is very large with respect to the 
temperature (i.e. the distribution is approximately not 
truncated), this new model yields the same results as the 
previous release and recapture model. 

The experimental results are shown in Figure [51 We 
repeated the experiment several times for different min- 
imum trap depths U m - ln . The data were recorded over 
a span of several days, where all other parameters were 
kept approximately constant. For each minimum trap 
depth we measured the recapture probability, P rr (Ai = 
0), the temperature, T, of the truncated Boltzmann dis- 
tribution and we calculated the associated mean energy, 
(E), of this distribution. The small day to day variations 
in the experimental parameters result in a small spread 
of the temperatures for a given value of C/ m i n . Never- 
theless, Figure [5jc) shows that the temperature of the 
Boltzmann distribution remains approximately constant 
down to very small trap depths, even when the distribu- 
tion becomes truncated. 

In contrast, the mean energy is reduced markedly as 
the distribution becomes more truncated, as we will now 
describe. The mean energy of a truncated Boltzmann 
distribution in a three dimensional harmonic potential 
can be expressed as: 

ftW E 3 e -E/k B T dE 
jUimnc E2 e -E/k B T dE 

' 1- (l + r] + W + hAe-v' 

= 3fc B T K —, i-J (3) 

1 - fl + ?7+ §r7 2 Je-*? 

where rj = Ut^ac/k^T . In the limit where the Boltzmann 
distribution is approximately non-truncated, i.e. rj 1, 
the mean energy is given by (E) — Sk^T. However, as 
shown in Eq. ([3]) the mean energy is reduced significantly 
as the truncation is increased, i.e. 77 — + 0. The probability 
that the atom survives this truncation process is given by 
Eq. @ as -P S urv(^trunc)- Naturally, this probability drops 
quickly as r\ — > 0, thus reducing the efficiency (or effective 
duty cycle) of this scheme. 

Figure [D(d) shows the mean energy calculated using 
Eq. §3§ with the measured values of the temperature and 



the truncation energy C/trunc- The solid line is a theoreti- 
cal prediction using Eq. ([3]) for a temperature T = 33 fiK. 

For minimum trap depths above 10 fiK, i.e. when the 
energy distribution is approximately non-truncated, the 
mean energy is (E) ~ 100 fiK. For minimum trap depths 
below 10 [iK, we see that the mean energy is significantly 
reduced. For example, for a minimum trap depth C/ m i n ~ 
0.4 fiK, we measure a mean energy of the atom at the 
end of this filtering sequence of (E) ~ 30 /xK, where the 
final trap depth is Uf ~ 2.8 mK. The truncation of the 
Boltzmann distribution thus leads to a reduction of the 
mean energy by a factor of 100/iK / 30/xK ~ 3.3. At 
the same time, the probability of recapture is reduced to 
P„(At = 0) - 0.1. 

In this section we have shown that if one is willing to 
compromise on the probability of recapturing the atom 
(or correspondingly, the duty cycle of the experiment), 
truncating the Boltzmann distribution is a good way to 
reduce the mean energy of the atom. 



VII. COOLING A SINGLE ATOM BY 
ADIABATICALLY LOWERING THE TRAPPING 
POTENTIAL 

Let us consider an atom in a trapping potential. By 
adiabatically lowering the trap depth the occupation 
probabilities of the vibrational levels are preserved [28| . 
For our optical dipole trap this conservation of the vi- 
brational number results in T/VU being a constant. In 
this section we will experimentally demonstrate that by 
slowly lowering the trap depth, the temperature of the 
atoms follows this relation, thus reinforcing our assump- 
tions about adiabaticity. 

To test the adiabaticity of our lowering scheme we per- 
form the same experimental sequence as in section IVI| 
but we now implement the "release and recapture" se- 
quence when the trap is shallowest, i.e. we measure the 
temperature of the atoms at the minimum trap depth 
f/min [ijl] . This is in contrast to section IVll where the 
temperature of the atoms was measured after the trap 
depth was returned to its initial value. 

A schematic of the experimental sequence is shown 
in Figure O together with the recapture probability, 
P„(At = 0), and the temperature, T, measured for vari- 
ous minimum trap depths U m i n . The data in figure [6](c), 
shows evidence that the temperature of the Boltzmann 
distribution follows closely an adiabatic lowering behav- 
ior. For comparison, we have plotted the adiabatic scal- 
ing law T/y/U = const, where the constant is equal to 
Ti/VUi with Ti = 33 fiK and Ui = 2.8 mK. 

The criteria for adiabaticity states that the rate of 
change of the oscillation frequency Co must be lower than 
uj 2 (uj is the oscillation frequency) at all times during 
the evolution [25|. Care must be taken to fulfill the 
adiabaticity criteria when the trapping potential is very 
small, corresponding to a small oscillation frequency. The 
rate of change of the oscillation frequency must be suffi- 
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FIG. 6: (Color online) Investigation of the adiabaticity during 
the change of the optical trap depth, (a) Schematic of the ex- 
perimental sequence (see text), (b) Probability of recapture 
for zero release time P rr (At = 0), as a function of the mini- 
mum trap depth (7 m in. The solid line is a theoretical fit to the 
data. Inset shows a schematic of the temporal experimental 
sequence, (c) Temperature T of the Boltzmann distribution 
plotted as a function of U m i u . The solid line is a theoretical 
prediction of the adiabatic scaling law T/VU = const. The 
dashed line is the extension of the solid line to minimum trap 
depths at which the density of states deviates from that of a 
three dimensional harmonic trap. 



ciently small in this case. In our experiment, this is done 
by smoothly ramping up and down the intensity of the 
dipole trap laser with acousto-optic modulators. The be- 
havior of the measured temperature in Figure E{c) gives 
us a strong indication that indeed the adiabaticity crite- 
ria is fulfilled when the trapping depth is small. 

Figure [S^b) shows the probability of recapture for zero 



release time P rr (Ai = 0) plotted as a function of the 
minimum trap depth E/ m i n , together with a theoretical 
fit. We reemphasize that for minimum trap depths less 
than t/ m i n ~ 10 /iK, the Boltzmann distribution becomes 
truncated. This implies that we are no longer in the 
three-dimensional harmonic approximation, and hence 
the density of states is no longer described by p(E) cx E 2 . 
Nevertheless, for clarity in Figure we have presented 
the results assuming that p(E) cx E 2 for all data points. 

As a typical result, we measure a temperature of 
T = 1.69 ± 0.02 /iK for the atoms in a trap that has 
been adiabatically lowered to a depth of U m \ n ~ 12 /iK. 
In this case, the Boltzmann distribution is only slightly 
truncated, implying that the quadratic density of states 
should still be valid. This temperature corresponds to 
a one hundred fold decrease in the temperature of the 
atoms compared to those that were directly loaded from 
the molasses. 



VIII. DISCUSSION 

We have shown that the single atom cooling via filter- 
ing techniques is quite efficient. Starting from an initial 
mean energy of 100 /iK after the laser cooling phase, we 
get a final mean energy of the atom of 30 /iK with a prob- 
ability of ~ 0.1. On our experiment, the loading rate of 
single atoms is on the order of 1 per second. If we are 
willing to comprise our duty cycle, we could obtain on 
average one atom every 10 seconds with a mean energy 
of 30 /iK in a ~ 2.8 mK deep trap, which corresponds to 
a ratio of the trap depth to mean energy of ~ 93. 

Another parameter characterizing the external motion 
of the atom in the trap is the mean vibrational number 
along one axis, defined as the ratio of the mean energy of 
the atom along this axis to the energy difference between 
vibrational states. As the trap has two different oscil- 
lation frequencies along the axial and radial directions, 
this corresponds to two different mean vibrational num- 
bers along these directions. If we assume the equiparti- 
tion of the mean energy among the three axis, we get for 
the mean vibrational number along the radial direction 
n± = (E)/3hv±. For the mean energy of 100 /iK that 
we obtained after laser-cooling the atom, this leads to 
iix ~ 4.3. By contrast, the lowest mean energy of 30 /iK 
that we obtained after filtering leads to n±_ ~ 1.3, which 
is 3.3 times smaller than the mean vibrational number 
before filtering. This mean number, close to 1 along the 
radial direction, raises some questions about the validity 
of our classical treatment for the temperature measure- 
ment for the lowest trap depths. This analysis, albeit 
interesting, is beyond the scope of this paper. This also 
places us in a good position to further cool the atom 
down to the ground state, for example by using Raman 
sideband cooling 

In the framework of quantum information, there has 
been proposed a protocol to entangle two atoms based 
on the emission of a single photon by one of the two 
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trapped atoms |29| . A crucial feature of this proposal is 
the ability to localize each atom within a distance smaller 
than the wavelength of the emitted photon, the so-called 
Lamb-Dicke regime. Reference [29| shows that for atoms 
described by a thermal distribution at temperature T, 
this condition becomes 77th = t?ld V 'k^T '/ 'huj <C 1, 
with to being the mean oscillation frequency of the atom 
in the trap, and 77LD = \J E r /huj being the Lamb-Dicke 
parameter (E r is the recoil energy). In our case, for a trap 
depth of ~ 2.8 mK, we find ?7ld = 0.20 and r] th = 0.56 
for a temperature of 33 /xK, as measured above. We are 
therefore in the right regime to apply this entanglement 
protocol. We finally note that 77th can be lowered when 
increasing the trap depth adiabatically. In this case ijth 
is related to the trap depth U by 77th °c U^ 1 ^ 4 . 

IX. CONCLUSION 

As a summary of this paper, we have investigated the 
energy of a single atom trapped in a tight optical dipolc 
trap. We have found that after applying laser cooling of 
the atom in the tweezer, the energy of the trapped atom 
follows a Boltzmann distribution. We have extracted the 



temperature of the atoms from the measured energy dis- 
tribution. We have also described a "release and recap- 
ture" method to measure the temperature of the atoms 
assuming a truncated Boltzmann distribution. This sec- 
ond method is in good agreement with the first one. 

We have shown that by adiabatically lowering the trap 
depth we have reduced the temperature of the atoms 
down to several microKelvin. Finally, we have cooled 
the atom by filtering the hottest atoms over many real- 
izations of the same experiment, the equivalent of evapo- 
rative cooling but at the single atom level. This leads to 
an atom close to the ground state in the radial direction 
of the trap. 
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